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Abstract 

A new probabilistic representation is presented for solutions of the incompressible Navier- 
Stokes equations in with given forcing and initial velocity. This representation expresses 
solutions as scaled conditional expectations of functionals of a Markov process indexed by the 
nodes of a binary tree. It gives existence and uniqueness of weak solutions for all time under 
relatively simple conditions on the forcing and initial data. These conditions involve comparison 
of the forcing and initial data with majorizing kernels. 



1 Introduction. 



The 3-dimensional Navier-Stokes equations governing the velocity of incompressible fluids with 
viscosity u and initial velocity uq are given by 

du ^ . ^ 

— — h u • Vu = vl\u — \/p + g 

ot 

V-ii = (1) 

uo{x) = u{x,0) 

Here u : x [0, oo) is the velocity vector, p : R^ x [0, oo) ^ R is the pressure and g : R^ x 

[0, oo) R^ is the external forcing. Although these equations have been studied extensively over 
the last century important open questions remain concerning existence and uniqueness of smooth 
solutions. Extensive background and development of results can be found in Temam ^H], Doering 
and Gibbon [Hj, Foias, Manley, Rosa and Temam 0, and Lemarie-Rieusset Jl]. Ladyzhenskaya 
jl2j provides a recent interesting and accessible review of results. 

*AMS 1991 subject classifications. Primary 35C15,60H30; Secondary 60J85. Key words and phrases: tree-indexed 
Markov process, branching random walk, incompressible Navier-Stokes equations. 
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There are two classes of existence and uniqueness results for solutions to the Navier-Stokes 
equations; existence and uniqueness for ah time if the forcing and initial data are small enough in 
some specified function space, or existence and uniqueness for a short time for larger initial data. 
This paper falls into the former camp in giving a probabilistic representation in physical space of 
the solutions of when a smallness condition is satisfied by the given forcing and initial velocity. 
This results in a method of demonstrating existence and uniqueness of weak solutions (in the sense 
of Lemarie-Rieusset |14j ) for all time under straightforward bounds on forcing and initial data. The 
focus is strictly on the Navier-Stokes equations in R'^. Indeed the method developed relies heavily 
on an integral formulation of the Navier-Stokes equations in R'^ incorporating Leray's projection 
onto the subspace of divergence free vector fields. 



The stochastic representation given here has parallels with the representations given by McKean 
jl5j for solutions of the Kolmogorov-Petrovskii-Piskunov equation and Le Jan and Sznitman 
for solutions of the Fourier transformed Navier-Stokes equations. In the seminal paper McKean 
gives the solution of the KPP equation 

du 1 d'^u 2 
with initial data u{x,0) = f{x), < f{x) < 1 for all x G R, as 

N{t) 
i=l 

where W = {W^'^ : i > 1} is a branching Brownian motion in R with initial value x. Roughly 
here one can think of the Laplacian term ^ as producing the Brownian motion, the term 
yielding the binary branching, and the —u term resulting in the exponential waiting time between 
branches. The existence of the solution is a consequence of the bound imposed on the initial data. 

LeJan and Sznitman in give the solution of the Fourier transformed Navier-Stokes equations 
in R3 

n(e,t) = e--l«l'*no(e) (2) 
+ ^ (2^)-3/2e-H«P-| _ j ^ u{7], t- s)^-ui^-rj,t- s)di] + ¥fg{i, t - 

in the form of a scaled expectation 



if both 



l"»(f)l and (3) 



are uniformly small enough. Here M is a multiplicative functional of a stochastic cascade rooted 
at ^ at time t with inputs derived from initial data and forcing scaled as in @ above. The 
scalar function H : R^ — > (0, oo] is given by either H{^) = or H{0 = a|C|-^e"°l«l for 

some a > 0. The vector projection onto the plane perpendicular to ^ has eliminated the 
pressure p. The term u{7], t — s)(, ■ u((^ — rj,t — s) yields binary branching in space according to the 
normalized Markovian kernel the coefficients e~'^l^l^* result in exponential waiting 
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times between branches. This general framework is elaborated and extended via the introduction 
of Fourier multiplier majorizing kernels in Bhattacharya et al j2| and [HI • Again the existence of 
solutions is a consequence of the functional bounds imposed on the forcing and initial data. 

The goal of this article is to construct a binary branching process with jumps that corresponds 
naturally to a formulation of solutions in physical space to the Navier-Stokes equations. It is not too 
surprising that this stochastic representation is more complicated than representations of solutions 
to the KPP equation or the Fourier transformed Navier-Stokes equations. Majorizing kernel pairs 
(/i, h) are introduced to be used as scaling multipliers. As such they both dominate the forcing 
and velocity and allow the construction of transition densities. Additionally, combined with the 
Laplacian term vAu, their presence can give rise to a /i-Brownian motion. The binary nature of the 
branching derives from the term u ■ Vn. In this representation, the waiting time between branches 
is related to first passage times for Brownian motion rather than being exponential in nature. This 
somewhat complicated relationship is in part due to the incompressibility constraint. This leads 
to the construction of a Markov process indexed by the binary tree rather than a Markov process 
in time. Analogous to the work of LeJan and Sznitman and Bhattacharya et al [2], the 
existence of solutions is a consequence of the bounds imposed on the forcing and initial data via 
majorizing kernel domination. The following theorem is representative. 

Theorem 1.1 Let h : — > (0,cx)] and h : — > [0, oo] with h locally square integrable and h, h 
jointly satisfying 

/ h^ix - y)\y\~'^dy < h{x) and / h{x - y)\y\^^ dy < h{x) for all x G R^. (4) 

If for all X G R'^ and t > 0, 

{4TTuty^/^\ [ uo{x-y)e~\y\^'^''*dy\<'Kvh{x)/ll and \g{x,t)\ < {'Kvfh{x) /II, 

then there exists a collection of probability measures {Px ■ x € R^} defined on a common measurable 
space {Q, J-), a measurable function H : (0, oo) xQ ^ R'^, and a weak solution u{x, t) to the Navier- 
Stokes equations (CJ) with 

Px{{lu ■.\E{t,uj)\ <2ttu /II for all t > 0}) = 1 for all x € R^ 

and 

u{x,t)=h{x) / E{t,uj)dPxiuj) for all xGR^,t>0. 
Jn 

Furthermore, the solution u is unique in the class {v £ {S'{R^ x (0,00)))^ : \v{x,t)\ < 2T:vh{x)/ll 
for all X £ R^,t > 0}. 

The probabilistic representation gives the scaled solution of the Navier-Stokes equations as the 
expectation of : t > 0}, which is defined in Sections 4 and 5 as a function of a Markov process 
indexed by the nodes of a binary tree. The scaling, h, is a non-negative majorizing kernel that con- 
trols the magnitude of the initial data uq. The kernel h is paired with a second kernel h that controls 
the magnitude of the forcing g. The conditions stipulated by (jlj) allow the definition of stochastic 
transition probabilities for an underlying branching Markov random walk while controlling the mag- 
nitude of the function H. For example, two majorizing kernel pairs satisfying are (/i(x), h{x)) = 
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(7r-3|x|-\(27r2|x|(l + |x|))-2) and {h{x),h{x)) = {Tr-^{l + \x\)-^,{4TT'^\x\{l + \x\)^)-^). This repre- 
sentation also gives a connection between Brownian motion in and solutions to the Navier-Stokes 
equations. 

The results of this paper are largely self-contained. The organization is as follows. Section 2 gives 
a definition of the solution space involved along with the derivation of an integral formulation of the 
incompressible Navier-Stokes equations incorporating Leray's projection. The integral formulation 
is given explicitly in terms of the transition density of Brownian motion in R'^. Majorizing kernel 
pairs are defined and some important properties of classes of majorizing kernels are derived in 
Section 3. Section 4 gives a probabilistic representation of the scaled weak solution to in the 
case of excessive majorizing kernels. This representation involves Doob's /i-Brownian motion in 
R^. This section also contains existence and uniqueness results in the case that the initial data 
uq is dominated appropriately by an excessive kernel h. An alternate probabilistic representation 
for general classes of majorizing kernels is given in Section 5. The existence and uniqueness results 
given in this section are consequently more general as well. The proof of Theorem 11.11 is included 
in Section 5. Section 6 contains a few concluding remarks. 



2 An integral reformulation of the incompressible Navier-Stokes 
equations. 

The delineation of the definition of solutions and the function spaces in which solutions exist is as 
follows. Throughout elements of R^ are treated as 3 x 1 vectors. The usual vector dot product and 
Euclidean norm | • | are used. ^(R^) = Cg°(R^) and ^(R^ x (0,cx))) = C^{B? x (0,oo)) where 
is the space of functions with compact support. 5'(R^ x (0,oo)) is the space of tempered 
distributions on R^ x (0, oo). 

Definition 2.1 A weak solution of the Navier-Stokes equations with initial velocity u{x,0) = uo(x) 
is a vector field u{x,t) G (5'(R^ x (0,oo)))'^ satisfying the following: 

(i.) u is locally square integrable on R^ x (0,cxd), 

(a.) V • n = 0, and 

(Hi.) there exists p £ S'(R? x (0, oo)) with ^ + u ■ Vn = z^An — Vp + g. 

After incorporating incompressibility via the Leray projection V and applying Duhamel's prin- 
ciple, the pressure term can be dropped and the Navier-Stokes equations written in integral form 
as 

u = e^^^UQ- f e'^'^^-'^^VV ■{u®u){s)ds+ f e''^^-'^^Vg{s)ds (5) 
Jo Jo 

V -uo = 

where the Leray projection V onto the space of divergence free vector fields of a vector- valued 
function v is defined as 
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See Foias et al jH| or Lemarie-Rieusset |14| for discussion. The integral formulation above is ex- 
pressed explicitly in Proposition 12. II that follows. 

Definition 2.2 A function h : ^ R is uniformly locally square integrable if for all G D(R^), 

sup / \(p{x - y)h{y)\'^dy < oo. 

Definition 2.3 A weak solution u of the Navier-Stokes equations on R'^ x (0, oo) is uniformly 
locally square integrable if for all (p G V^R? x (0, cxd)) 

sup / / \(p{x — y,t) u{y,t)\'^dydt < oo. 

xeR3 JR3 J(0,oo) 

The following equivalence theorem can be found in Lemarie-Rieusset |14| . 

Theorem 2.1 If u is uniformly locally square integrable, then u is a weak solution of (0j if and 
only if u is a solution of 



The following definitions are useful in stating Propostion l2.1l Let 

K{y, t) = (27rt)~3/2g-|s,p/2t^ ^ > ^ ^ 1^3 

denote the transition density of a 3-dimensional Brownian motion starting at the origin at time 
t = 0. For y G R^ with \y\ > 0, let Cy = y/\y\ denote the unit vector in the direction of y and let 
Fy denote the 3 by 3 matrix that projects vectors onto the space perpendicular to y; so that the 
entries of Py are 

Define the bilinear forms bi and b2 on R^ \ {(0,0,0)} x (R^)^ via 

bi (y; u, v) = {u- ey)'PyV + {v ■ ey)'PyU 

and 

h2{y] n, v) = bi(y; u,v) +u - {I - Se^^e* )i; Cy. 



Proposition 2.1 If u : R^ x (0,oo) — > R^ is locally square integrable and satisfies 



u{x,t) = / uq{x - y)K{y,2vt)dy 

JR3 



+ /" ( \ \--^K{z,2us) hi{z-,u{x — z,t — s)^u{x — z,t — s)) (6) 

+ f -j— 2z^s) — . f M / K{y,2us)dy'] h2{z;u{x — z,t — s),u{x — z,t — s)) 

^'^\^\ J{y-\y\<\z\} ' 

^-p7(/ - Se^e*) / K{y,2vs)dy\g{x - z,t- s)\dzds, 

W hy\y\<\A^ ' J 



+ ^/^(.,2..)P.-^^| 
then u is a weak solution to 0). 
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This integral formulation incorporating incompressibility is key to the proof of Theorem 11.11 
and the other existence and uniqueness results given in Sections 4 and 5. The following lemmas 
are used in the proof of the proposition. Throughout the Fourier transform on R'^ is defined as 



/(O = (2vr)-3/2 [ e-"<f{x)dx 

J-R? 

and understood in a distributional sense. Then also 

f{x) = (2vr)-3/2 f e--«/(e)de. 

J-R? 

The bilinear forms bi and b2 defined above and the matrix / — Se^e* can be bounded as follows. 



Lemma 2.1 For any y G with \y\ / and u,v e R^, 

\hi{y;u,v)\ < \u\ \v\, 
\h2{y;u,v)\ < 2\u\ \v\, 



and 



\{I-3eyel)u\ < 2\u\. 



(7) 
(8) 

(9) 



Proof: Fix y,u,v € R^ with \y\ / 0. Let a = Cu ■ Cy, f3 = Cy ■ Cy, and 7 = e„ • e^. For a, 6 G [0, 1], 
/(a, b) := a(l — 6^)^/^ + 6(1 — a^)^/^ < 1 with equality for + 6^ = 1. The triangle inequality gives 

\hi{y;u,v)\ = \u\v\ {aPye^ + /3Pj/e„| < \u\ \v\f{\a\, |/3|) < |u| 

To derive (jSJ, first notice that 

|bi(y; = |M|>p(a2(l - /J^) + 2ap{-f - a/3) + /^^(l - a^)) 

and 



Then 



|7 - aP\ = \eu ■ Pj,e„| < |Pj,e,| |Pe,| = (1 - a^f'^{l - P^f'^. 



\h2{y;u,v)\^ = \hi{y;u,v)\'^ + {u ■ {I - ?,eyei)vf 



< 



< 



u\' \v\' {a' + + a'P' + (7 - 2a(3y) 

n|2 \v\^ (a2 +p2^ ^2^2 ^ ((^ _ ^2)l/2(^ _ ^2)1/2 ^ |^^|)2) 

n|2 (1 + a^p^ + 2|a/3|/(|a|, (1 - P^)^/^)) 
\v\^ {1 + \ap\f 



< {2\u\\v\f. 

To verify @ note that 

I (/ - 3eyel)u\^ = u-{I + 3eye\)u = \u\^{l + Sa^) < 4|iip. 



□ 
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Lemma 2.2 Suppose that for each fixed s > 0, T{-, s) : R^^^ has Fourier transform 

t{^,s) = (27r)-3/2e-H?P^p^. 

Then for x E R^ with \x\ ^ 0, 

T{x,s) = K{x,2iys)P^ - {4tt)~^\x\'^{I - 3e^ei) / K{y,2us)dy. 

J{y-\y\<\A} 

The result presented in Lemma [2.2l is not new; it goes back at least to Solonnikov ^7j, who used it to 
calculate estimates of solutions of the linearized Navier-Stokes equations. It appears again in Koch 
and Solonnikov [Tl] in a study of the Stokes problem. More recently, Thomann and Guenther ^2] 
use it to compute an explicit formula for the fundamental solution of the linearized Navier-Stokes 
equation in terms of special functions. The representation given here is more explicitly stated in 
terms of the transition density for Brownian motion, than that of Solonnikov |17j . Koch and 
Solonnikov or Thomann and Guenther |19| . The derivation is included for completeness. 

Proof: Fix s > and notice that k{i, s) = (27r)-3/2g-i.|a2,^ Suppose that 70(6 = (27r)-3/2|^|-2e-''l«l- 
Some computation gives 

7o(x) = (27r)-3/2(2j,s)-i/2|^|-i e-^'/^^'dp. (10) 

Jp=0 

(For computational details of the Fourier transform for radial functions, see Folland [Hj, page 247.) 
If 7jj' : R^ ^ R is defined via 

= -(2v^)-3/2pe-HeP^ 

then 

= (27r)-3/2(2z.s)-i/2|<5„-(|xr2e-l-lV4-_|^|-3 e-pV4..^^) 



p=0 



+XiXj\x\-^ (3\x\~^ e-p'/^'^'dp - (31x1-2 ^ (2i.s)-i)e-l^l'/4'^^) | . (11) 
Combining (fTU]) and (|TT|) gives 



Tij{x,s) = K{x,2i^s)6ij + jij{x) 
= K{x,2iys){P^)ij 



-(27r)-3/2(2z.s)-i/2|^|-3 (■ y e-^'/^'-'dp - \x\e-^^^" [6ij - 3x,Xj\xr' 

To complete the proof note that 

/ K{y,2us)dy = {TTiys)-^/^( e-p'^^'^'dp - \x\e-\''\'/^''') . (12) 

J{y:\y\<\x\} ^Jp=0 ' 

□ 

The following lemma is given without proof. 
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Lemma 2.3 Suppose that V : R — > R has Fourier transform 



Then V is given by 



V{C) = i{2Tiy^'^ / v{r]) C-v{^- 7])dr]. 



3 Q 

Vkix) = —{Vk{x)Vj{x)). 
j=l J 



Lemma 2.4 For V as defined in Lemma W^ 

I ^iy,s)V{x - y)dy = - \ {4:i^s)~'^\y\K{y,2i's)hi{y;v{x - y),v{x - y)) 

■/r3 Jr3 

+ (\y\-'K{y, 2us) - 3(47r)-i / K{z, 2us)dz) b2(y; v{x - y),v{x - y))}dy. 

^ J {z:\z\<\y\} ' ^ 

Proof: Lemma 12.21 gives 

^i,j{y,s) = Sijgi{\y\) +yiyjg2{\y\) 

where for r > 

gi{r) = (47ri/s)-3/2e-'^V4^« _ (4^)-V-3 f K{y, 2vs)dy 

J{y-\y\<r} 



and 



c,2(r) = 3(47r)-V^^ / K{y, 2us)dy - {Attvs) 

J \v-\y\<r} 



'{y-\y\<r} 

A term by term integration by parts leads to the vector-valued function 

T{y,s)V{x -y)dy = / \{2usy^ K{y,2us)[{ey ■ v{x - y)f y - y ■ v{x - y) v{x - y)] 

+ 92{\y\){^{y ■ v{x - y))v{x -y) + {\v{x - y)\^ - h{ey ■ v{x - y)f)y^ jrfy 
= / \ -{^i'sy^\y\K{y,2us)hi{y\v{x-y),v{x-y)) 
+\y\92{\y\)b2{y;v{x - y),v{x - y))^dy. 

□ 

Proof of Proposition \2. 1{ The incorporation of Leray's projection into the Fourier transformed 
Navier-Stokes equations as seen in gives a convenient formulation for weak solutions n of 
For fixed s, apply the inverse Fourier transform to the Fourier Navier-Stokes equations Q, letting 
V represent the vector-valued function with Vfc(x, s) = X]j=i a§"(^fc(^> s)uj{x, s)), and use Lemmas 
TTR ESI and El to see that 

u{x,t) = / uo{x-y)K{y,2vt)dy 

J-R? 
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+ 



/ {~ '^iy,s)V{x - y,t - s)dy + T{y,s)g{x - y,t - s)dy\ds 
Js=o [ Jus Jr3 J 



uo{x - y)K{y,2vt)dy 

t 



+ j yj 3 y-^'^^'> \y\I^{y^ 2i/s) bi(y; u{x -y,t- s),u{x -y,t- s)) 

+ (\y\~^K{y,2iys) -3{iTr)-^\y\-^ [ K{z,2vs)dz] 
^ J{^--\A<\y\} ' 

b2(y; u{x - y,t - s),u{x - y,t - s)) ] dy 



+ J^^ r(y, s)g{x -y,t- s)dy^ds. 



□ 



3 Navier-Stokes majorizing kernel pairs. 



The majorizing kernel pairs defined below are used to control the magnitude of the forcing and the 
initial velocity in a way that allows the control of the magnitude of the velocity over time. The first 
component of the pair, h, dominates the velocity. The second component, h, dominates the forcing 
in a way tailored to maintain the velocity domination. The conditions specificing the majorizing 
kernels allow construction of transition probabilities based on the integral representation of the 
Navier-Stokes equations given in Proposition 12. II above. These kernels play a role in physical space 
analogous to the role played in Fourier space by the scaling functions of Le Jan and Sznitman jj^] 
and the Fourier multiplier majorizing kernels introduced by Bhattacharya et al in j2j. 

Definition 3.1 The pair {h, h) : x R'^ (0, oo] x [0, oo] is a majorizing kernel pair with constant 
pair (7,7) G (0, 00) x [0, 00) if h is lower semi- continuous and uniformly locally square integrable, 

/r3 h\x-y)\y\~^dy 

sup -Ei — — = 7 < 00, (13) 

xeR3 h{x) 

and 

/r3 Hx - y)\y\''^dy 

sup -E^^ — = 7 < CX). 

XGR3 Hx) 

Jf 'J = J = we say that (h, h) is a standard majorizing kernel pair. 

Notice that both h and h are finite a.e. with respect to Lebesgue measure and h[x) = a.e. if and 
only if 7 = 0. A simple scaling argument shows that if a majorizing kernel pair (/i, h) has constant 
pair (7,7) with 7 > 0, then {'-/~^h, (77)"^/i) is a standard majorizing kernel pair. 
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Recall from potential theory that an excessive function h : (0, oo] satisfies both 



/r3 h{y)K{x - y, 2vt)dy a v f 

sup — < 1 and lim / 

=R3,t>o Kx) no 7r 



h{y)K{x — y, 2ut)dy = h{x). 

xeR3,t>o riyx) *iu iR3 

See Bass 1^ or Doob X, fo'^ background. Majorizing kernel pairs of particular interest are those 
with the first member h being excessive. 



Definition 3.2 A majorizing kernel pair (h, h) is an excessive majorizing kernel pair if h is ex- 
cessive. 



The initial data uq and the forcing g will be compared, respectively to h and h. 



Definition 3.3 The pair {uo,g) = {uQ,{g{-,t) : t > 0}) is {h,h)- admissible if 



sup < oo ana sup — < oo. 

xeR3 H^) xeR3,i>o h{x) 



Example 1: Let ho{x) = \x\ ^ and let h{x) = hQ{\x\) with j-^^ h{x)dx = 47r J^Qr'^hQ{r)dr < oo. 
Then (/iq, h) is an excessive majorizing kernel pair with constant pair (vr^, 47r^ /r=o f'^hQ{r)dr). This 
can be seen by observing that /iq is uniformly locally square integrable and for x 7^ 0, 



/ \y\-^K{x-y,t)dy = \xr^ P{\Z\ <\x\/Vt), 



/R3 

where Z is a standard normal r.v., 

|y|^^|x — yl^'^dy = 7r'^\x\^^ , 

/R3 

and ^ 

\x\ / \y\~^h{x — y)dy = Air / min{r, |x|} r ^o('^)'^^- 

JR3 Jr=0 

The value on the left above increases to Air f^Qr'^hQ{r)dr = f-^a ho{\y\)dy as \x\ — > 00. As an 
example, let ho{r) = r~^(l + r)~(^+^) for some fixed e > 0. Then {uo,g) is (/iq , /i)-admissible if 
suPx.gR3 klho(a^)| < 00 and sup^.gR,3^t>o + i)| < 00. 

The Fourier transform of the function /iq given here is the Le Jan-Sznitman scaling function 
1^1 ~^ that arises naturally in the study of the Fourier Navier-Stokes equations. See Le Jan and 
Sznitman for the original study. See also Bhattacharya et al 2 for refinement of the Fourier 
space majorizing kernel method. See also the remark following Example 4. 

The following propositions illustrate some of the structure of the class of majorizing kernel pairs. 
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Proposition 3.1 Let {h,h) be a standard (excessive) majorizing kernel pair. 

(a) Then 

{h{- - fi),h(- - fi)), /xGR, 
{ah{a-),a^h{a-)), a > 0, 

and 

{h{A-), h{A-)), A a 3 by 3 matrix with A^A = I, 

are also standard (excessive) majorizing kernel pairs. 

(b) Let F denote a probability distribution function on R^. Then 

if h{--y)dF{y)J h{--y)dF{y)) 
is also a (excessive) majorizing kernel pair with constant pair (7,7) G (0, 1] x (0, 1]. 

Proof: Part (a) is easily checked via the appropriate change of variables. Part (b) follows from an 
application of the Cauchy-Schwartz inequality followed by Fubini's Theorem. □ 

Proposition 3.2 Let {{hj,hj) : j > 1} be a sequence of (excessive) majorizing kernel pairs with 
corresponding constant pairs {(7j,7j) '■ j > !}• 

(a) Then {hi A h2,hi A /i2) is a (excessive) majorizing kernel pair with constant pair (7,7) G 
(0,71 A 72] X [0,71 A 72]. 

(b) For any p G (0, 1) 

{hlhl-^,hlhl-n 

is a (excessive) majorizing kernel pair with constant pair (7,7) G (0,7i72~^] x [0,72*72"^]. 

(c) For {pj : J > 1} with "^"^iPj = 1 and pj > for each j > 1 

00 00 
j=i 3=1 

is a (excessive) majorizing kernel pair with constant pair (7,7) G (0, ^°^]^Pj7j] x [0, ^°^;^Pj7j]. 

Proof: Part (a) is obvious. Apply Holder's inequality to deduce (b). Part (c) again follows from 
the Cauchy-Schwartz inequality and Fubini's Theorem. □ 

Definition 3.4 The functions hj : —>■ (0, 00], j = 1,2 are equivalent if for some c G (l,oo), 

c~^/ti(x) < h2{x) < chiix) for all x G R^. 

Proposition 3.3 If hi and hi are equivalent and {hi,h) is a majorizing kernel pair, then {h2,h) 
is also a majorizing kernel pair. 
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Proof: The proof follows easily from the definitions. 



□ 



Example 2: The majorizing kernel /iq of Example 1 has a singularity at the origin. Propositions 
3.1 and 3.2, can be used to construct a kernel with a countable number of singularities as follows. 
Take {fij : j > 1} to be a sequence in and {pj : j > 1} with Y^^=iP3 — ^ ^'^d Pj — each 
j > 1. Let 

oo 

h{x) = ^Pj\x - Hjl'-^ 
i=i 

and 

oo 

H^) = ^Pjho{\x - 

for some fiQ : [0,oo) [0,oo] with J^Qr'^ho{r)dr < oo. Note that if fiQ has a singularity, then h 
itself will have a countable number of singularities as well, {h, h) is a majorizing kernel pair with 
constant pair (7,7) G (0,vr^] x (0, /j^3 /lod^^Ddx]. 

Example 3: The operations of Propositions 3.1 and 3.2 can also be used to produce bounded 
kernels based on /iq. For example, convolving /iq with the probability density f{y) = (27r|y|)~^(l + 
yields the bounded kernel 

hi{x) = (i + ixr^. 

Taking h{x) = J-^^ ho{\x — y\)f{y)dy yields a majorizing kernel pair with constant pair (7,7) S 
(0,7r3] X (0,/j^3 hoi\x\)dx]. 

Using Proposition 13.31 we see that any radial function h that decreases as |x| increases with 
lim|2.|_^o = c G [0,oo) and lim^^^^^ \x\h{x) = C S (0,oo) can be used as a majorizing kernel 

when paired with an appropriate h. Indeed kernels h of this description are in the Marcinkiewicz 
space see Cannone and Karch ,5,^ for a related study of stability results for Navier-Stokes 

equations. 

Example 4: Let H{x) = (1 + |xp)^^. Although H is not excessive, it is uniformly locally square 
integrable and 

/ H{y)K{x - y, t)dy < (1 + ?,e-^^'^)H{x) for ah x G R^. 

Additionally, 

/ H'^{y)\x -y\"^dy = TT^H{x) for ah x G R^ 
However, it is not difficult to show that for any h with Jj^3 h{y)dy > 0, 

svip{H{x))^^ / \y\^^h{x — y)dy = 00. 
xeR» ^R3 

This gives the majorizing kernel pair {H,0) with constant pair (7r^,0). Thus {uo,g) being {H,0)- 
admissible corresponds to sup^.gj^3(l + |j;|^)|Mo(a^)| < 00 and the forcing g being identically 0. Note 
that H{(,) = \/ 7r/2|^|~-^e~'^L This is the second of the two Fourier Navier-Stokes majorizing kernels 
introduced by Le Jan and Sznitman [l3] . 



12 



Remark: Although Examples 1 and 4 are suggestive, it is not generally true that the Fourier 
transform of a majorizing kernel h satisfying (|13|) is a Fourier multiplier majorizing kernel or vice 
versa. The key to the correspondence in Examples 1 and 4 is the equality 



for a fixed constant c. The condition h * /i(^) < C\(,\h{(,) is that required by Bhattacharya et al [2] 
for Fourier multiplier majorizing kernels for the Fourier transformed Navier-Stokes equations. 

Example 5: For p G (1,2], let Hp[x) = (1 + Then {Hp,0) is a majorizing kernel pair with 

constant pair (7, 0), 7 < 7r^~'P(l + l/\/2)P~^. This can be deduced by first noting that H as defined 
in Example 4 is equivalent to (1 + with c = 1 + l/\/2 and then applying (b) of Proposition 

1,3.21 using H and the kernel (1 + of Example 3. 



4 Construction of the stochastic representation and existence and 
uniqueness with excessive kernels. 

This section opens by introducing the notation necessary to define tree-indexed Markov processes. 
After defining some crucial transition densities, the integral equation (jH)) of Proposition 12.11 is 
reformulated in terms of the transition densities and the velocity and forcing scaled respectively 
by h and h for a majorizing kernel pair {h,h). The underlying tree-indexed Markov process and 
stochastic recursion are then defined and described. The existence and uniqueness results for 
excessive kernels follow. 

Let V := U^o{0, 1}" denote the fuh binary tree with root = {0,1}°. Let dV := {0,1}^ 
denote the boundary of V. For v =< vi, ...,Vn >G {0, 1}", we say that the magnitude of v, \v\ = n. 
For f G V define u|0 = (p and for v ^ (j) and n > 1, v\n =< fi, >. If for some n < \v\, 

v\n = w, we say that w is an ancestor of v and, conversely, that t; is a descendant of w. For 

V G V with = n, < n < 00, set v = v\{\v\ — 1) and v * k =< vi, k >, k = 0,1. Thus 

V is the child of v and has children v * and v * 1. Although it may currently seem nonsensical, 
it is convenient in the sequel to interpret (p as the precursor of the root (j). For v,w £ V with 

V ^ w, denote the last common ancestor of v and t;; as w A tn. This is defined as follows. Let 
Uy^ii, = max{m > : v\m = wlm} and then take v Aw = We say that W C V is a rooted 
binary sub-tree if € W, for any v S W, v\k £ W for all k < \v\, and for any v * j £ W, j = 0, 1, 
then V * {1 — j) € W as well. For a finite rooted binary sub-tree W, define the boundary of W, 
dW, as the elements of W which have no descendants in W; that is 



The interior of W, W° = W \ dW, consists of the elements of W that have descendants in W. For 
example, a Galton- Watson tree with a single progenitor and offspring distribution concentrated on 
and 2 is a rooted binary sub-tree. If the expected number of offspring is less than or equal to 1, 
then the binary sub-tree is finite with probability 1. 




On the Fourier side, this becomes 



dW = {v eW ■.v*O^W}. 
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Let (B,;B) be a measurable space and let X = {Xy : v E V} he a V-indexed collection of B- 
valued random variables defined on a common probability space. Let J^^ = a{Xy), = 0"({X^|„ : 
n < \v\}) denote the cr-field generated by the collection of r.v.'s indexed by v and ?;'s ancestors, 

and 7iy = a{{Xu, : w \v\ = v}) denote the cr-field generated by the collection of r.v.'s indexed by v 
and v's descendants. Let J^j denote the trivial cr-field. 



Definition 4.1 'X. is a V-indexed Markov process if for any v E V with \v\ < oo, 0.^*0 and 7iv*\ 
are conditionally independent given T-o, so that 

P{AinA2\J'y) = P{Ai\J^y)P{A2\Ty) a.s. P for Ak e ny,k,k = 0,1, 

and for any v,w eV and any T-L^ -measurable random variable Y with E\Y\ < oo, 

E{Y\g^) = E{Y\J^^^„) a.s. P. 



The distribution of a V-indexed Markov process is completely specified by the conditional distri- 
butions of the Xy given (ot v eV. 

The integral form of the incompressible Navier-Stokes equation will be expressed in terms of 
an expectation of a functional of a V-indexed Markov process. The following conditional densities 
involving the majorizing kernel pairs will be used to specify the transition probabilities and thus 
the distribution. Let {h, h) be a majorizing kernel pair. Fix x G and define a probability density 
on X R^ via 

^ \ynz\-^h\x-z)m>\y\] 
J[y,z\x) 27r J^,\z\-''h^x-z)dz ' ^ ' 

If h is not identically 0, likewise define 

ZTT Jj^3 \z\ ^ri[x — zjdz 

The above densities will provide spatial transition densities for our V-indexed Markov process. 
Define two conditional waiting time densities for s > and y G R^ via 

fo{s\y) = (27r)-V2(2^)-3/2,-5/2|y|3g-|,|V4.. = 27rs-'\yfK{y, 2us) (16) 

and 

fi{s\y) = (47ri.)-V2s-3/2|y|e-l^lV4.^« = ^7ru\y\K{y, 2us). (17) 
Adjoin a trap state to R^, and define 



' ih{x))-^h{y)K{x-y,2i^t) if x,y £ , h{x) < oo; 

J(y t\x) = \ ^~ (^(^))^'/r3 h{y)K{x - y, 2vt)dy if xeR?,y = 9; 

|1 \i X = y = 6 

otherwise. 

If h is excessive, this is the transition density of a /i-Brownian motion on R^ U {6}. 



(18) 
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For a majorizing kernel pair (h, h) and x S U {^}, define 

Xoix) = I Q ^ ^ ^ (19) 

and 

^(rr,t) = 4^. (20) 
h{x) 

Also let 

, ^ /r3 \y\-'^h\x-y)dy 
mix) = — , (21) 



and 



_ /R3|y| ^Mx-y)rfy . . 



Proposition 4.1 Suppose that the majorizing kernel pair {h,h) is used to define f, f, /o, /i, J, 
Xo; V'; CLnd m as given in ^14] ) through h2l3\) . If h{x)x{x, t) is locally square integrable and x 
satisfies 



X{x,t) = / X(i{y)J{y,t\x)dy (23) 

+ / \m{x)[ fo{s\z)f{y,z\x)hi{z;xix - z,t- s),xix - z,t- s)) 

Js=o Jj/eRS J^eR3 I V 

+ (2fi{s\z)f{y,z\x) - 3fi{s\y)f{y,z\x)^h2{z;xix - z,t- s),x{x - z,t- s)) 
+m{x) (2fi{s\z)f{y, z\x)Pz(pix - z,t- s) 

-h{s\y)hy,z\x){I - 3e^e*)(/?(x - z,t- s)) \dzdyds 



then u{x,t) = h{x)x{x,t) is a weak solution to 0). 



Proof: Suppose that x{x,t) satisfies above and define u{x,t) = h{x)x{x,t). Note that 

|z|^^/i^(x — z) 



and 



Then 



f{y,z\x)dy = 
yGR3 jjis\y\ ^h\x-y)dy 

7, I V, \z\-^h{x-z) 

j[y,z\xjay = — = . 

yeR3 j-^3\y\~^h{x -y)dy 



u{x,t) = / uo{y)K{x - y,2ut)dy 
JyeR^ 



f f I hix)m{x) ( /o(s|2;) 

+ / / St — I I 2^.2/ \T\ 112 ^\KZ\u{x - z,t- s),u{x - z,t- s)) 

Js=oJz&iv^ [ \y\ ^ - y)dy \ kr 
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+ f ~ / ^T^irTd^y) ^2(^;; u{x - z,t- s),u{x - z,t- s)) 

Jr3 |y| - v)dy V fI 



\y\<\z\ 2vr|y||2; 
/ uo{y)K{x-y,2ut)dy 



s=0 JzeR' 



+ 11 {[ -r^ hi{z;u{x - z,t - s),u{x - z,t - s)) 



Aus 



+ _ ^ / 2us)dy) h2{z; u{x - z,t - s),u{x - z,t - s)) 

\ \z\ 47r|z|* J|j;|<|^| / 

^^'(z, 2z^s) P^fi'(x — z,t — s) 
The result then follows from Proposition 12.11 □ 



It is helpful to combine the spatial and temporal transition densities /, /, /o and /i given in 
equations H14|) through ()17() via the following randomization. Let {pk : 1 < k < 5} satisfy pk > 
and Ylk=iPk = 1- To appropriately balance the randomization, fixp G (0, 1/2] and take pi = p/ll, 
6p/ll, and P4 = Ps = (1 —p)/2. Note that Pi + P2 + P3 = P ^ 1/2. For fixed 



Pk 

P2 = 4p/ll, P3 
xGK^ let 



f{y,z,s,k\x) 



PiMs\z)f{y,z\x) if k = l 

P2fiis\z)f{y,z\x) if A; = 2 

P3fiis\y)l{y,z\x) if A; = 3 

P4/i(sk)/(y,2|x) if A; = 4 

P5/i('S|2/)/(y,^|2;) if A; = 5. 



(24) 



This can be thought of as the joint conditional density of a quadruple {Y, Z, r, «;) taking values in 
R3 X R3 X (0,oo) X {1,2,3,4,5}. 



We now define the tree-indexed Markov process that underlies our probabilistic representation. 
Let X = {({KW : t > 0},X„y, ,Z^,Ty,Ky) : V G V} be a V-indexed Markov process with the 
following transition probabilities. First assume that, given that Xjj = x, {Vv{t) : t > 0} is a h- 
Brownian motion with initial value x and transition density J as defined in p8j) above. Note that 
each Vy is itself a Markov process which can be assumed to be continuous with probability 1; see 
Doob 0, Part 2, Chapter X for background. Take the transition density of (Yy , Zy , , k^) given 
that Xy = X to be given by f{y, z, s, k\x) as defined in (j^U, and then set 

Xy Xy Zy. 

Finally assume that for all f € V, the process {Vy{t) : t > 0} and the ensemble {Xy,Yy, Zy,Ty, Ky) 
are conditionally independent given Xy. 
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Notice that the distribution of : t ^ 0}, X^, Yy, Z^, r^, only depends on the ensemble 

{{Vjj{t) : t > 0},Xjj,Yjj, Zfj,Tjj, Hjj) through the r.v. Xjj. The ensemble {X^ : f € V} is itself 
a branching Markov random walk. The distribution of X then depends on the initial value X^. 
Denote the probability measure corresponding to = x by Px and the expectation with respect 
to this probability measure by E^. 

For V & V, define the random functional as follows. Let denote the random bilinear 
operator 

, j bi(Z,„; • , • ) if = 1 

■ ' ■ ^ " I {-ir^h2{Z,; • , • )/2 if K, = 2,3 ^^^^ 

and denote the random matrix 



_j Pz, if K„ = 4 



Then set 



^vit) =Xo{Vv{t)) + ^^""^^^^ B^(T^,o(t - r^),T^.i(t - r^))l[/^^ = 1,2,3] n [r^ < t] 

+ CMXv,t - r,)lK = 4, 5] n [r, < t]. (27) 

Notice that the recursion implicit in the definition of terminates if Ty > t or «;„ = 4 or 5. That 
is 

T,(t)l[r, >t]= Xo{Vv{t)) (28) 

and 

T^{t)l[K^ = 4, 5]1[t^ <t]= Xo{Vv{t)) + ^^^^ CMXv,t - r^). (29) 

The random space-time branching mechanism giving rise to can be explained as follows. 
Initialize the Markov process X by taking X^ = x G R"^. Activate the /i-Brownian motion V0(i) 
starting at at time t = 0. This motion will run until time t with Xo(^(i)) being used in 
calculating the first term of T^plt). If has entered the trap state 9 this term has value 0. 
Although the process runs until time t, if is less than t the path is deactivated at time r^. If 
deactivation occurs before time t, there are two possibilities. If /t^ = 4 or 5, the forcing is input into 
T<^ by evaluating C^(/? at location and time t — t^. If = 1, 2 or 3, two new active paths Vq 
and Vi are started at location X^ at time r^. The distributions of these newly activated paths are 
conditionally independent given X^. If the paths Vq and Vi have been activated, then the process 
repeats on each of them with X^ taking the place of and t — taking the place of t. This is 
illustrated in Figure 1. Notice that all paths activated by time t are used in calculating T0(t), even 
though they may be deactivated before time t. For example evaluating the stochastic recursion for 
Figure 1 gives 

T<^(t) = Xomt)) + H!^^B^(To(t - r^), Ti(i - r^)) 

with 

To(t - r^) = Xo(Fo(t - r^)) + ^^^Co<^(Xo, t-r^- tq) 
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Figure 1: An illustration of an initial segment of the random process X with k,^ = 1, 2 or 3, kq = 4 
or 5, and ki = 1,2 or 3. The deactivated paths are represented by dotted lines. 



and 

- r^) = Xo(Vi(i - r^)) + H^^^^Bi(xo(V^i,o(i - - n)),xo{Vi,i{t - - n))). 
Here t is evaluated at the end of the time interval illustrated. 

Remark: Although it is not obvious from the construction above, the waiting times can be 
tied to first passage times for Brownian motion. Here's a sketch of an alternate construction that 
highlights the relationship. For simplicity, just consider a single f G V. Let : t > 0} be a 

Brownian motion in with W(0) = and transition density K{y,2vt), independent of and 
the ensemble {{V^{t) : t > Q],X^,Y^,Z^,k^). Write W{t) = {W^^\t),W^^\t),W^^\t)) and for 
a > let 

r(^)(a) = inf{s>0:W^W(s) = a} for z = 1,2,3. 

The r'^*)(a)'s are independent stopping times with individual densities {^ttu)"^/'^ s~^/'^ae~"^ ^^^^ . 
Let ^ 

T^'\a) = (£{T^^{a)r'r\ 
1=1 

It is easy to check that r(°)(a) has density (27r)-V2(2i/)-3/2s-5/2^3g-aV4i^s. Note that T^^) is not 
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a stopping time with respect to the usual filtration associated with the Brownian motion W. Let 

r r(o)(|z,|) if = 1 

= <^ r«(|Z,|) if = 2,4 

[ r«(|y,|) if = 3,5. 

It is easy to check that, conditioned on the value of X;^, the joint distribution of {{V^{t) : t > 
0}, Xy,Yy, Z^,fj;, K^), IS the same as that of : t > 0}, Xy,Yy, Z^,Ty, k^)- This allows us to 

think of the r^j's as waiting times that reflect excursion levels in orthogonal directions of a hidden 
independent Brownian motion in R^. 



Theorem 4.1 Let {h, h) he an excessive majorizing kernel pair. If T<^ is as defined in \2'1\) and 

Ex\T^{t)\ < M < oo for all x and t, then 

u{x,t) = h{x)ExT^{t) 
is a weak solution to the Navier-Stokes equations. 

Proof: To(t — r,^) and Ti(t — t^) are conditionally independent given J~(j) , giving 
EA'^4t)\^^) = Xo{V4t)) 

+rn{x)!^p^^hi{Z^-,Ex{To{t-T^)\J^^),Ex{ri{t-T^)\J^^))l[K^ = 1] 

+ (2p^H[K^ = 2] - 3ip3y'l[K^ = si 

b2(Z^; Ex{To{t - T^)\T^),E,{Ti{t - t^)\T^)) ^1[t^ < t] 

+m(x)|2p4 1i[k^ = A]Pz, - {p^r'lK = 5](/ - 3ez,e*^^)|^(x - Z^,t - t^)1[t^ < t]. 

Since also 

Ec,{Ti{t - T^)\r^)l[Z^ = z,T^ = s\= E^-.T^t - s) for i = 0,l, (30) 

E^T^it) = Ex{E,{T4t)\T^)) 
J{y,i\x)xo{y)dy 



R3 



+ 



/ / / <m{x)[fo{s\z)f{y,z\x)hi{z;Ex^:,T^{t-s),Erc-zT^{t-s)) 

Js=0 JyeR3 JzeRS [ V 

'^fiis\z)f{y,z\x) - 3fi{s\y)f{y,z\x)^h2{z;Er,-zT^{t - s),Ex_^T^{t - s))^ 



+ m(x)( 2fi{s\z)f{y,z\x)Pz(p{x - z,t- s) 

fi{s\y)f{y,z\x){I - 3ezel)(p{x - z,t- s)] )dzdyds. 



Then u(x,t) is uniformly locally square integrable since = h{x)\ExT^{t)\ < Mh{x). From 

Proposition 14. II the function h{x)ExT^{t) is a weak solution to (0). □ 
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Theorem 4.2 (Existence and Uniqueness) Let {h,h) be an excessive majorizing kernel pair with 
constant pair (7,7) G (0, S-TTi/p/ll] x [0,27rz^(l — p)] for some p £ (0,1/2]. If {uQ,g) is {h,h)- 
admissihle with 

l«o(a:)| . , \9{x,t)\ 
sup < ae and sup — < pe 

x&R? H^) xGR3,t>0 h{x) 

for some a G [0,1) and e, P £ (0, 1 — a), then u{x,t) = h{x)ExT ^{t) is a weak solution to the 
Navier-Stokes equations with 

\u{x,t)\ 

Furthermore, this solution is unique in the class {v G (5'(R^ x (0, 00)))^ : sup^.gj^3 f>o ^^y^^)^ ^ 

The conditions of Theorem 14.21 are sometimes more easily checked in the fohowing form. 

Corollary 4.3 Let {h,h) be an excessive majorizing kernel pair with constant pair (7,7), 7 > 0, 
and suppose that {uQ,g) is (h, h)-admissible with 

\uo{x)\ STTupae \g{x,t)\ {4TTjyfp{l-p)(3e 
sup , , , < — — and sup < — — r 

X6R3 h{x) II7 x&RM>0 h{x) II77 

for some p G [0, 1/2], a G [0, 1) and e, /? G (0, 1 — a), then u{x, t) = h{x)ExT^{t) is a weak solution 
to the Navier-Stokes equations with 

\u(x,t)\ S-Kupe 

This solution is unique in the class {v G (5'(R^ x (0, 00)))^ : supj.gR3,t>o ^^h(x)^ — ^ ii^'^ I • 

Proof of Corollary 14-^ ' Rescale {h,h) and consider the majorizing kernel pair 

8^ (4vr.)Ml-p) ^) 
4l7 ' II77 ^ 

with constant pair (Svrz^p/ll, 27rz^(l — p))- □ 
The next corollary treats the case of no forcing. 

Corollary 4.4 Let (h, 0) be an excessive majorizing kernel pair with constant pair (7, 0) and sup- 
pose that 

\uo(x)\ Airvae 

for some a G [0,1) and e G (0, 1 — a), then u{x,t) = h{x)ExT fj,{t) is a weak solution to the 
Navier-Stokes equations with 

\u{x,t)\ AiTve 
sup , , . < ——. 

xGR3,t>0 Kx) II7 

This solution is unique in the class {v G (5'(R'^ x (0, 00)))'^ : sup^gR3^t>o ^^h(x)^ — Tl^i- 
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The proof of Theorem 14.21 depends upon Proposition 14.21 that foUows. 



Proposition 4.2 Let {h, h) be an excessive majorizing kernel pair with constant pair (7, 7) € 
(0, 8iTiyr]p/ll] X [0, 27ri^ry(l — p)] for some p G (0, 1/2] and 77 > 0. Suppose that for some a G [0, 1) 
and e, 13 G (0,(1 — a)/r]), {uo,g) is [h, h) -admissible with 

\uo{x)\ , , \9{x,t)\ 

sup , , < ae ana sup — < pe. 

xeR3 n[x) xeR^t>o h{x) 

Then for all x G and t > 

|T<^(t)| < e a.s. Px. 
The following lemma is used in the proof of the proposition. 

Lemma 4.1 Let W C V be a finite binary sub-tree. Suppose that {hy : x R'^ —>■ R'^ : v G W} 

has 

sup |b^(j;, y)| < |x| |y| for all x,?/GR^, 
and {{Uv, ZviVvjCTv) G R^ x R^ x [0, cxd) x {0, 1} : ?; G W} satisfies 



and 



sup rjy < T] < 00 



sup \yv\ < ae and sup |z^| < /3e for some a G [0, 1) and e,/3 G (0, (1 — a)/r]]. 



Then x^ defined iteratively on W via 



X7- 



satisfies 



Uv + o-yTjyZy if V e dW 

yv +rjvKixv*o,Xv*i) if v£W°. 



sup < e. 



Proof: The result follows by induction starting on dW. If u G dW, then 

\xv\ < \yv\ + i]v\zv\ < ae + rjPe < e. 
If u G W° and \xy^k\ < e for k = 0, 1, then 

I Xy I ^ I yv I ~i~ 1 1)^ {xy^o, Xy^i)\ < ae + r/e^ < e. 

□ 

Proof of Proposition Let W = {v £ V : Ky^ = 1,2,3 for all j < \v\}. Then W is a random 
binary sub-tree. Indeed, W corresponds to a Gallon- Watson tree with each individual having either 
or 2 offspring. The probability of 2 offspring is p = P{Ky = 1, 2, 3) < 1/2 and the probability of 
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offspring is 1 —p = P{nu = 4,5). With probability 1, W is finite; see for example Harris iTIT. Set 

= X G and = 0. For i; G W with \v\ > 0, let = Efclo^T^lfc- For fixed t < oo, T^{t) 
is a functional of the random ensembles indexed by the nodes of the a.s. finite binary sub-tree 
>V(t) C W defined by 

W{t) = {v eW : <t}. (31) 
Recalhng (jSHJ and ((221) , if u G dW{t), 

T^{t - S^) = xo{V^{t - S.o)) + ^'^^_^^^ CMXv,t- - t^)1[k^ = 4,5] n [r^ < t - 
From (EZI), for v G W°{t) 

+ ^ — ^B^(T^^,o(t — Sv — Ty),Ty:ti{t — Sv — Ty))l[Kv = 1, 2, 3]l[r^ <t — S^]. 

P 

This is the setting of Lemma HH] with y„ = XoiVvit — S^)), = B^, Zy = Cyip{Xy,t 



( if vedW{t) 

if v(^w°{t) 



and ay = l{Ky = 4,5]. Lemma [2.11 gives |B^(j;,y)| < \x\\y\ for all x,y £ R^. By assumption we have 
IxoCKjI-s))! < ae for all s, and again using Lemma 1^11 \Cyip{Xy, s)\ < fie for all s. The multipliers 
rjy are bounded as follows: 

llm(x) II J^i\y\-'^h'^{x -y)dy II7 

= — < < rj 

p S-Kuphyx) Sttup 

and 

4m{x) /r3 \y\~^h{x - y)dy , 1 , 

— ^ — = ^ ^ s 

1— p ^-KviX — p)h{x) 2iTi'{l — p) 

for all x G R^. Lemma l4. II gives |T(^(t)| < e a.s. P^- From Theorem 14. II then 

\u{x,t)\ = h{x)\E^T^{t)\ < h{x)E^\T^{t)\ < eh{x). 

□ 

Proof of Theorem \4-.'fl\ The existence of the weak solution u{x,t) = h{x)ExT ^[t) follows immedi- 
ately from Theorem 14. II and Proposition 14.21 with ij = 1. Uniqueness is derived using a martingale 
argument as follows. Suppose that u is a solution to © with 

\v{x,t)\ 
sup < e. 

xeR3,t>o '^[^J 

Set p{x,t) = v{x,t)/h{x). Take = x G R'^ fixed, fix t > and let 

>V(")(t) = G W(t) : |v| < n} (32) 
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where W(i) is as defined in ()31|). For n > and v G define the random functionals ^1"^ via 

the iterative construction 

= XoiKit)) 

+ ^ ^""^^^ ^ B^{p{X^,t-T^),p{X^,t-T^))l[K^ = 1,2,3] n [T^<t] 

+ ^""^^^^ C^y?(X^, t - r^)lK = 4, 5] n [r^ < t] 
1 — p 

and, for n > 1, 

= Xo(K(t)) 

+ ii^^B4M/,!:-^)(t - r,),^ti'\t - r.))l[K. = 1,2,3] n [r. < t] 
4mfX-) 

+ —-^^Cv^{X,a,t-T^)l[Ky=4,5]n[Tv<t]. 
1 — p 

Notice that for each n, ^'^"^(i) depends only on those ensembles in X indexed by f G W^"^(t). If 
n;(")(*) = then ^'J'^(t) = T^{t). Since is a solution to ®, 

= h{x)p{x,t) = h{x)E.Ml^\t) 

and for all 

(t - s,)\g^)i[t - 5, > 0] = - - 5, > 0]. 

If for some n > 1 and all ?; G V, ""^^(t — Sy)\Qv) = p{X^,t — Sy) on the set [t — S,; > 0], then 

- s,)\g,)i[t -s,>o] = xoiVvit - s,)) 

+ {{nm{X^)/p)B,{E,{¥Zo'\t -Sv- T,)\g,),E,{^^Zi'\t -Sv- T.)ia.))l[K. = 1,2,3] 
+ {Am{X^)/{l - p))CMXv,t - - t^)1[k^ = A,5]}i[t^ < t - S^] 

+ ^{llm{Xij)/p)Bv{p{Xv,t - Sv - Tv), p{Xv,t - Sv - Tv))1[k.v = 1,2,3] 
+ {Am{X^)/{l - p))CMXv,t - - t^)1[k^ = A,5]}i[t^ < t - S^]. 
and, on the set [t — 5^ > 0], 

J{y,t- Sy\Xjj)xo{y)dy 



/R3 

+ / / / \m{X^)lfo{s\z)f{y,z\Xij)hi{z;p{Xy-z,t-Sy-s),p{Xy-z,t-Sv-s))) 
Js=o Jr3 Jr3 I V 

+ i'^fiis\z)fiy,z\Xij) - 3fi{s\y)f{y,z\Xii))h2iz;p{Xij - z,t - Sy - s),p{X^ - z,t - - s)) 



+ m{Xv)[2fi{s\z)f{y,z\X^)PMXv -z,t-S,-s) 

-fi{s\y)f{y,z\Xy){I - 3ezel)(f{Xy - z,t- Sy- s)^ ^dzdyds 

= p{Xy,t-Sy). 
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By induction then, Ex{"^v {t — Sy)\Qjj)l[t — > 0] = p{X^, t — Si,)l[t — 5^, > 0] for all n and v. In 
particular, Ex'^lp\t) = p{x,t) for all n. Since yV{t) is finite a.s. Px, there exists with probability 1 

some random with W^"'\t) = W{t) and = T(^(t) for all n> N. The random functionals 

and are all bounded by e in magnitude, so 

\p{x,t) - x{x,t)\ = \Ex^l;;\t) - ExT4t)\ < Ex\^'-^\t) - T^(t)| < 2eP,(w(")(t) / W(t)). 
This probability goes to as n goes to oo. □ 

Example 1: Take p = a = l/2to see that if 

sup |x||no(x)| < iz/llvr^ and su]i\g{x,t)\ < hQ{\x\) for all x G 

for some ho : [0, oo) [0, oo] with J-^^ ho{\x\)dx < /IItt, then there exists a unique weak solution 
u{x,t) to (P) with sup^, J i)| < 2^'/ll7r^. 

Example 2: Again take p = a = l/2. If there exists {fij : j > 1} in and {pj : j > 1}, pj > 0, 
= 1 with 

|uo(x)| < (z^/llvr^) ^^Pj|x — and sup \g{x, t)\ <^'^^pjho{\x — fij\) for all x G R^ 

for some ho with ho{\x\)dx < v'^/IItt, then there exists a unique weak solution u{x,t) to 
with supj |u(x,i)| < (2z^/ll7r^) Y^Pj\x — fJ-j\~^ for all x G R^. 

Example 3: If 

sup (1 + |x|)|uo(x)| < z^/llvr^ and sup |(7(x, t)| < /i(x) for all x G R^ 

for some h with < h{x) < J-^^ hoi\x — y|)(27r|y|)^^(l + \y\)^'^dy where ho{\x\)dx < v'^/Utt, 
then there exists a unique weak solution u{x,t) to with sup^.^t(l + |x|)|ti(x, i)| < 2z^/ll7r^. 

Notice that Examples 1,2 and 3 all give existence and uniqueness in a sub-space of the Marcinkiewicz 
space L^'°°; c.f. Cannone and Karch 0. 

5 Non-excessive kernels and an alternate probabilistic represen- 
tation. 

The following representation can be used with both non-excessive and excessive majorizing kernels. 
Let {h, h) be a majorizing kernel pair and let X^^^ = Yy, r^,, k^) : u G V} be the V-indexed 

Markov process with the transition density of (1^, Z^,, r„, k^,), given that = x, as defined in (|24|) 
and again let Xy = Xy — Z„. To define the stochastic recursion, let 

, ^\ /r3 ^o(a; - y)K{y, 2vt)dy 
mo{x,t) = -£i — 
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and let xo, V-, ""t-j and in be as defined in equations ((TO)) through ^F^. We will also use the random 
operators B^, and as defined in (|25|) and (|26j) . Now define the random recursive functional H 
on V via 

= mo(X{;,t) + ^^'^^'^'"^ 'By{~.^^o{t - Ty),E^^i{t - r^))l[Kt, = 1,2,3] n [t^ < t] 
P 

+ , ^ C^^{X^, t - T^)l[Ky = 4, 5] n [r^ < t]. 
1 — p 

If (/i, h) is excessive, this random recursion essentially replaces the first term of T^, Xo(^v{t))y with 
mQ{Xy,t), the conditional expectation of Xoi^vit)) given Qj). 

Proposition 5.1 Let (h, h) he a majorizing kernel pair with constant pair (7,7) G (0, STivrjp/ll] x 
[0,27rz^7/(l — p)] for some p E (0,1/2] and r] > 0. Suppose that for some a £ [0,1) and £,/? E 
(0, (1 — a)/r]), {uo,g) is {h, h)-admissible with 

\ j^suoix -y)K{y,2ut)dy\ \g{x,t)\ 

sup — — < ae and sup — < pe. 

xgR3,t>o K^) xGR3,t>o h{x) 

Then for all x E R'^ and t > 

|S</,(i)| < e a.s. Px- 

Proof: The proof of this proposition exactly parallels that of Proposition 14.21 Here the y^ of 
Lemma im is equal to mQ{Xy^t — Sy), so by assumption, \y^\ < ae. The proof is then identical. □ 

The conditions appearing in the following corollary are somewhat simpler. 

Corollary 5.1 Let {h,h) be a majorizing kernel pair with constant pair (7,7) E (0, Svrz/Typ/ll] x 
[0, 2ttl'7]{1 — p)] for some p £ {0,1 /2] and rj > and 

J^3h{x -y)K{y,2ut)dy 
sup — < M < 00. 

xeR3,t>o "-(^J 

Suppose that for some a £ [0,1) and e, P £ {0,(1 — a) /rj), {uo,g) is {h,h)- admissible with 

l^o(a;)| . , \9{x,t)\ 

sup < ae/M and sup — < pe. 

xgR3 h^x) xeR3,t>o h{x) 

Then for all x £ Y{? and t > 

\'^4>{t)\ < e as. Px- 

The proofs of the following theorems are omitted due to their similarity to proofs appearing in 
the previous section. 

Theorem 5.2 Let (/i, h) be a majorizing kernel pair. If Ex\'E^{t)\ < M < 00 for all x and t, then 

u{x,t) = h{x)E^E^{t) 
is a weak solution to the Navier-Stokes equations. 
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Theorem 5.3 (Existence and Uniqueness) Let {h, h) be a majorizing kernel pair with constant pair 
(7,7) G (0, Svrz^p/ll] X [0, 27rz^(l — p)] for some p £ (0,1/2]. If {uQ,g) is {h^h)- admissible with 

\ j^3UQ{x -y)K{y,2vt)dy\ \g{x,t)\ 
sup — < ae and sup — < pe. (33) 

a;GR3,t>0 "(2^) xeR3,i>0 h{x) 

for some a G [0, 1) and some e, /? G (0, 1 — a), then u{x,t) = h{x)Ex^if,{t) is a weak solution to 
the Navier-Stokes equations with 

\u{x,t)\ 
sup —rrv~ - ^■ 

This solution is unique in the class {v G (5'(R^ x (0, 00)))^ : sup2,gR^3^t>o ^^i^^^)^ ^ 



Three useful coroharies fohow. The first imposes separate bounds on the magnitude of the 
initial velocity relative to h and the size of h* K relative to h. In the last two corollaries a single 
bound is imposed on the magnitude oi uq * h relative to h. The last corollary treats the case of 
no forcing. This is of particular interest in regards to regularity; see remark (4) in the following 
section. 



Corollary 5.4 Let {h, h) be a majorizing kernel pair with 

/r3 h[x - y)K{y,2ut)dy ^ ^ 

sup — — < M < 00. (34) 

xeR^,t>o "-(^J 

and constant pair (7,7) G (0,87rz^p/ll] x [0, 27ri/(l — p)] for some p G (0, 1/2]. // {uo,g) is {h,h)- 
admissible with 

\Mx)\ , . , \9ix,t)\ 
sup < ae/M ana sup — < pe (35) 

xGR^ n(x) a;GR3,t>0 h{x) 

for some a G [0, 1) and e, /3 G (0, 1 — a). Then u{x,t) = h{x)Ex'^^{t) is a weak solution to the 
Navier-Stokes equations with 

\u{x,t)\ 

This solution is unique in the class {v G (5'(R^ x (0, 00)))^ : sup^gj^3^^>o ^""^^-^^ < ^l- 

Proof of Corollary \5.4[ Combining (|34|) and (|35|) gives (|33|) . □ 
The next corollary relies on the rescaling seen in the derivation of Corollary 14.31 

Corollary 5.5 Let {h,h) be a majorizing kernel pair with constant pair (7,7), 7 > 0, and suppose 
that (no, (7) is {h,h)- admissible with 

\ !^iUQ{x - y)K{y,2vt)dy\ Snupae \9{x,t)\ ^ {4TTufp{l-p)f3e 
sup — < — — and sup — < — — : 

xGR3,t>0 II7 xGR3,t>0 h{x) II77 

(36) 
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for some p £ [0, 1/2], a £ [0, 1) and e, /? S (0, 1 — a), then u(x, t) = h{x)Ex'^^{t) is a weak solution 
to the Navier-Stokes equations with 

\u{x,t)\ STTupe 

xeR3,t>o Hx) II7 

This solution is unique in the class {v G (5'(R^ x (0, 00)))^ : sup^gj^3^t>o ^^h{x)^ — TT^i- 

Proof of Theorem Equation Q) guarantees that is satisfied with 7 = 7 = 1 and p = a = 

1/2. Then (jSZj) holds for any e < 1/2. □ 



Corollary 5.6 Let {h, 0) be a majorizing kernel pair with constant pair (7, 0) and suppose that 

I /r3 ^0(2; - y)K{y, 2vt)dy\ Anvae 

for some a € [0,1) and e G (0, 1 — a). Then u{x,t) = h{x)Ex'^(f,{t) is a weak solution to the 
Navier-Stokes equations with 

This solution is unique in the class {v G (5'(R^ x (0, 00)))'^ : sup^gj^3^t>o ^^h(x)^ — Tl^}' 



Example 4: If sup2,(l + |xp)|no(x)| < ll7r/16i^(l + 3e ^Z'^) and the forcing (7 is identically 0, then 
there exists a unique weak solution to (P) with sup^ j(l + |xp)|u(x,i)| < llir/Su. 

Example 5: Taking Hp{x) = (1 + \x\)~p for p G (1, 2] fixed, 

/ Hp{x - y)K{y, 2vt)dy < (1 + 2,e-^iy\l + l/^y/^Hp{x) 

JR3 

for all X G and t > 0. Then, if sup^(l + IxD^liiolx)] < 7rP-3z^(l + 3e-2/3)-p/2(i + i/^)i-3p/2/ii 
and the forcing g is identically 0, there exists a unique weak solution to (0) with s\iY>xti)- + 
\x\y\u{x,t)\ < 27rP-3i/(l + l/V^)i-P/ll. 

6 Remarks. 



(1) It appears that the integral formulation on which the probabilistic representation is based can 
also be used as the basis for Picard iteration. The advantage of the probabilistic approach as it 
appears here is that it allows, with probability 1, the assumption that there are a finite number of 
iterates. This permits easy derivation of bounds via inductive arguments. 

(2) The probabilistic representations of the Fourier transformed Navier-Stokes equations given 
by Bhattacharya et al i2j give both short time existence and uniqueness of solutions for large initial 
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data and decay over time of the Fourier transformed solution. It is possible that a physical space 
representation related to the ones given in this paper will also indicate short time existence and 
uniqueness as well as the time decay of solutions. 



(3) A semi-group approach underlies this probabilistic representation. Montgomery-Smith ^ 
uses such an approach to demonstrate blow-up in finite time of an equation related to the Navier- 
Stokes equations. A more complete understanding of the integration of the bilinear forms bi and 
b2 may aid in understanding possible blow-up of the Navier-Stokes equations in R^. 

(4) In some cases the results given in this paper can be used to demonstrate regularity and 
consequently uniqueness in a larger class of functions. The argument is as follows: assume that 

is regular and take the forcing to be identically 0. Take (/i, 0) to be a majorizing 
kernel pair with h G L2(R^) n Lg(R^) for some q > 3. If sup^ j | / uo{x — y)K{y,2vt)dy\/h{x) is 
small enough, then sup^ ^ \u{x-, t)\/h{x) < e for suitable e > and for any r > and T < oo 

j\ j \u{x,t)\UxYdt < ^'^'^ j\ j h'i{x)dxydt < oo. 

That is, the Ladyzhenskaya-Prodi-Serrin condition is satisfied, and the solution u is regular and 
unique in a larger function space; c.f. Temam |18j . Chapter III. A suitable kernel here \s h = Hp 
of Example 5 with p G (3/2, 2]. 
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